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SUMMARY 


Two rocket models having cruciform, aspect-ratio-5, unswept tails 
and a fineness-ratio-20 fuselage were flight tested over a Mach number 
r an ge of approximately 1.7 to 2.4. One of the models had cruciform, 
aspect-ratio-5 *4 forward surfaces in line with the tails. The models 
were given step disturbances by pulse rockets at intervals throughout 
the Mach number range and stability derivatives were obtained from the 
measured responses. The roll rates of the models varied from 10 radians 
per second to 5 radians per second. 

The measured lift-curve slopes in the plane in which the disturbance 
originated (pitch) were much lower than the lift -curve slope in the other 
plane (yaw) for both models. The measured lift-curve slopes in the pitch 
plane were in good agreement with potential-flow theory. The damping of 
the model with tail only was greater than that predicted by theory. The 
damping of the model with both forward surfaces and tails was about the 
same as that predicted by theory. The aerodynamic center of both models 
was farther forward than was predicted by potential -flow theory. 


INTRODUCTION 


As the speed of missiles has advanced far into the supersonic regime, 
the general tendency has become more and more to use wings of low aspect 
ratio. This trend is a logical one since the increase in lift-curve slope 
with aspect ratio has become very small at Mach numbers of 2 and above, 
and the elastic problems of a high-aspect-ratio wing at high speed become 
rather severe. As a result, very little experimental data has been 
obtained or needed on high-aspect-ratio wings in this speed range. 
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Recently, however, with the advent of folding-fin rockets as airplane 
armament, high-aspect-ratio wings at Mach numbers around 2 become of 
interest. The chord is restricted in order to allow the fins to fold 
readily and a comparatively large span is often necessary to have suf- 
ficient stabilizing area. It is conceivable that guidance systems may 
be developed which can be fitted in small aircraft rockets and will 
require a rather blunt nose body. Such a missile will probably have a 
high-fineness -ratio fuselage and forward-control fins since the rear part 
of the fuselage will be, of necessity, a rocket motor. 

Adequate theory exists for calculating the lift-curve slope and 
aerodynamic center of such a high-aspect-ratio configuration in a super- 
sonic potential flow. (For example, see refs. 1 and 2.) A question 
exists as to whether the assumption of potential flow is valid for 1 
conditions that may be encountered by guided aircraft rockets. The pur- 
pose of the present investigation is to compare stability derivatives 
obtained experimentally at high Reynolds number with stability deriva- 
tives calculated from existing theory for a possible folding-fin config- 
uration. Experimental static and dynamic stability data as obtained 
from two free-flight models are presented. The fuselages of both models 
had hemispherical noses and body fineness ratios of 20. The cruciform 
tail surfaces of both models were unswept and untapered and had an aspect 
ratio of 5. One of the models had cruciform lifting surfaces on the for- 
ward part of the missile such as might be used for control fins. The 
models were disturbed at predetermined intervals during their flights 
by small rockets firing normal to the model flight path. Dynamic and 
static stability derivatives as calculated from the measured responses 
of the missiles to these disturbances are compared with the derivatives 
as calculated from potential-flow theory. 


SYMBOIS 

A n normal acceleration, g units 

A^ transverse acceleration, g units 

C-p) drag coefficient, 

qjtr^ 

C m pitching moment about model center of gravity, moment 

1 qSc 

Cjy normal-force coefficient 

C T lift coefficient 
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resultant-force coefficient 


side-force coefficient 


Ixjlyjlz m0D ients of inertia about X-, Y-, and Z-axes, respectively, 
slug-ft 2 

% space motion factor (used in ref. 5) 

Pt total pressure, lb/sq in. 

R Reynolds number per foot 

V free-stream velocity, ft/sec 

X,Y,Z body coordinate axes 

( span of tail surfaces (used only in pb/2V), ft 
exponential damping constant in e"* 3 ^ per second 
c mean aerodynamic chord, ft 

g acceleration due to gravity, 32.2 ft /sec 2 

m model mass, slugs 

p roll rate, radians/sec 

q dynamic pressure, lb/sq ft 

S exposed tail area in one plane, sq ft 

t time from model launching, sec 

r body radius, ft 

x distance from center of pressure of air-flow indicator to model 

center of gravity (4.74 ft for both models) 

a angle of attack of model, deg 

ot^ angle of attack indicated by flow-direction indicator, deg 

P angle of sideslip of model, deg 

Pi angle of sideslip indicated by flow-direction indicator, deg 
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8 pitching velocity, radians/sec 

ijf yawing velocity, radians/sec 

oj damped natural frequency of resultant motion, radians/sec 

Subscripts: 

q derivative with respect to 0c/2V 

a derivative with respect to a 

a derivative with respect to acy 

P derivative with respect to p 

p derivative with respect to 

t indicates value at trim condition 


^Mrh) 

P 

^ c / 2V (j7^) 


A dot over a. or 0 indicates a derivative with respect to time. 


MODELS AND APPARATUS 


Model Description 

The models used for these tests had fuselage fineness ratios of 20 
and cruciform tails of aspect ratio 5 that were untapered and unswept. 

One model, referred to throughout this report as model 1, had no forward 
surfaces. The other model, referred to in this report as model 2, was 
identical in configuration to model 1 except for the addition of aspect- 
ratio-3.4 cruciform forward surfaces arranged in line with the tail. A 
general arrangement of both models is shown in figure 1 and details of 
the lifting surfaces are shown in figure 2. The telemeter antennas shown 
in figures 1 and 2 and the total -pressure tubes shown in figure 1 were 
indexed 45° with respect to the tails in an effort to minimize inter- 
ference. The air-flow direction pickups located in front of the models 
were believed to be located far enough forward to minimize their wind- 
shield effect on the blunt hemispherical noses. The fuselage forward 
section was made of 3/l6-inch wall-steel tubing and the rear section was 
a standard 5-inch HVAR rocket motor. The solid-steel fins were welded 
to the fuselage. It is believed that the models were sufficiently rigid; 
thus, effects of flexibility were negligible. Photographs of the models 
are presented in figure 3* The dimensional and mass data are as follows: 
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Exposed tail area (per plane), sq ft 0.6944 

Fuselage cross-sectional area, sq ft O.I 363 

Exposed forward wing area (model 2 only), per plane, sq ft . . 0.4l67 

Tail mean aerodynamic chord, ft 0.4l67 

Center-of -gravity position (both models), inches from nose . . 43.9 

Model 1 Model 2 

(fuel burned) (fuel burned) 


Weight, lb . 

.... 162.1 

168.6 

IX, slug-ft 2 

.... 0.293 

0.32 

Iy or Ig, slug-ft 2 

.... 38.9 

39.4 

£ 

iy 

.... 0.0075 

0.0081 


INSTRUMENTATION 


The models were equipped with an NACA eight-channel telemeter which 
transmitted a continuous record of normal, transverse, and longitudinal 
acceleration, pitch angular acceleration, roll rate, total pressure, angle 
of attack, and angle of sideslip. No reliable information was received 
from the pitch angular acceleration pickup on either model. 

Velocity was measured by a CW Doppler velocimeter and the positions 
of the models in space were measured with an NACA modified SCR 584 
tracking radar set. Atmospheric temperature, pressure, wind velocity, 
and wind direction were measured by a Rawinsonde. The method of refer- 
ence 3 was used in an attempt to launch the models on a day when the 
atmospheric turbulence was low. 


Test Technique 

The models were launched from a near-zero-length mobile launcher 
at an elevation angle of 45°. Each model was boosted to supersonic 
velocity by a first-stage booster powered by two 6-inch-diameter solid- 
propellant rocket motors which together produced approximately 12,000 
pounds of thrust for 3 seconds duration. The 5-inch HVAR rocket motor 
which made up the rearward part of each of the model fuselages was fired 
directly after first-stage burnout, separated the model and booster, 
and propelled the model up to a maximum speed. As the models coasted 
down through the Mach number range, four small rockets on each model were 
fired normal to the flight path at intervals in order to cause free pitch 
oscillations of the models. These small rockets, referred to as pulse 
rockets, delivered 20 pound-seconds of impulse over about 0.06 second. 


CONFIDENTIAL 



6 


* **C01^liEJnf!ML NACA RM L57A31a 

PRECISION OF DATA 


The velocity data, as obtained by the CW Doppler velocimeter, were 
corrected for flight-path curvature and for wind effects at altitude. 

The angle of attack a and angle of sideslip 3 at the model center 
of gravity were obtained from the wind-flow-direction indicator by using 
the following relationships: 



(37.3)(32.2)A n 

V 



3 = 


Pi + 


Pi 


- (57.3)(32.2)^ - 


Since all linear accelerometers obviously could not be located at the 
model center of gravity, corrections had to be applied for angular veloc- 
ities and accelerations. These corrections were very small, in almost 
all cases being less than 1 percent of the instrument range. 


ACCURACY 


The accuracy of the instrumentation used for these tests should 
cause the following quantities to be within the following incremental 
limits for the two Mach numbers listed: 


Mach number 

a - Oj. or 0 - 

Cjj or Cy 

C D 

p 

2.4 ± 0.02 

±0.15 

±0.003 

±0.05 

±0.2 

1.7 ± 0.02 

± 0.15 

± 0.005 

± 0.07 

±0.2 


The absolute values of a and 3 which are not used for results in this 
report are not nearly as accurate as a - and 3 - 3-^ because of 

possible fin misalinements of the air-flow indicator. No rigorous assess- 
ment of the accuracy of the stability derivatives can be obtained. In 
the method used for determining CN a > the random errors and some of the 

systematic errors are canceled out when the slope is determined. In the 
case of C^ and Cm^ + Cm^, the linearity of the derivative itself is 

often very questionable. The following listed accuracies are based on 
limited experience from instances when several similar models of a con- 
figuration have been tested or when checks with other methods have been 
possible: 
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Cj^ or C Y p, percent ±5 

Cq^, percent ±8 

Cmq + C^, percent ±20 


RESULTS AND DISCUSSION 


All coefficients presented in this report with the exception of Cp 
are based on the exposed area of the tail in one plane (0.6944 sq ft) 
and the tail chord (0.4l67 ft). All moment coefficients are referred 
to the model center of gravity (43*9 inches from model nose). The drag 
coefficient Cp is based on the fuselage cross-sectional area. The 
variation of Reynolds number per foot with Mach number for the tests is 
shown in figure 4. The measured total -pressure variation with time for 
model 1 is presented in figure 5 along with a calculated value obtained 
by using Rayleigh's supersonic pitot equation. The Mach number and static 
pressure used in the pitot equation were obtained from the CW Doppler an d 
SCR 584 radar and Rawinsonde measurements. When compared, the measured 
and calculated total pressures for model 2 appeared to be very much like 
those shown for model 1. For both models the agreement below a Mach 
number of 1.9 is very good but above a Mach number of 1.9 the measured 
total pressure is lower and very irregular. This result suggests that, 
above M = 1.9, the total -pressure tube is in a region of separated flow 
and it seems logical that at other positions around the fuselage there 
may also at times be a separated flow. Attempts to correlate the irregu- 
larities in the measured total pressure with variations in a and (3 
showed nothing that appeared systematic but they were evident even at 
angles of attack of less than 1°. Occurrences of this nature, of course, 
are strongly dependent on the Reynolds number of the test. 


Time Histories 

The -time history of %, Cy, roll rate, and Mach number during pulse- 
rocket firings is shown in figures 6 and 7 for models 1 and 2, respec- 
tively. Note the irregular response of C^ and Cy because of their 
being referenced to a conventional body-axis system while the model is 
rolling at from 5 to 10 radians per second. Little evidence of cross 
coupling is apparent on the roll-rate trace because of the small angles 
of attack and sideslip (usually less than 2°). The rather abrupt steps 
on the roll-rate trace occur when the pulse rocket fires, since the thrust 
line of the pulse rockets was not in the X,Z plane but was parallel 
to it and displaced 0.75 inch. The average variation of pb/2V with 
Mach number is shown in figure 8 to be approximately constant for model 2 
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but increasing with Mach number for model 1. The roll rate could not 
be accurately predicted from fin misalinement . From measured wing mis- 
alinements, model 1 was calculated to have had a value of pb/2V 
of -0.0021 during flight test. Similarly, the effects of wing-tail 
interference being neglected, model 2 was calculated to have had a value 
of pb/2V of 0.0010. 


Normal Force Due to Angle of Attack 


Plots of Cn against a are presented in figure 9( a ) for model 1 
and in figure 9(b) for model 2. Similar plots of Cy against 0 are 

presented in figure 10(a) for model 1 and in figure 10(b) for model 2. 

Summary plots showing Cm and Cy against Mach number are presented 

a p 

in figure ll(a) for model 1 and in figure ll(b) for model 2. Theoretical 

values for Cm or -Cy as calculated from reference 1 (solid line 

% 1 0 

curves) are also shown in figures 11(a) and ll(b). As may be seen in 
figure 11, the agreement between the experimental points of C^ and 

the theoretical curve is good for both models but the experimental points 
for -Cyp are considerably higher than the theoretical values for both 

models. Since this result did occur on both models it is very unlikely 
that it is due to any instrument errors. It is believed that, if for 
some reason the vertical tails were in a more turbulent flow than the 
horizontal tails, the experimental values of Cy^ would be less than 

theory rather than greater. It should be noted that the experimental 
values of C 


values of C 


(or Cy 'j are total derivatives, whereas the theoretical 


to account 


for 


a\ 

presented are partial derivatives. Attempts were made 
Cflj due to § and Cy due to ijr but no appreciable 


difference in the derivatives CL T and -C v resulted. A somewhat 

N a Y p 

similar phenomenon occurred during the test reported in reference k 
where C was the same for a model that rolled at about 5 radians 

per second as for a model that did not roll if the pitch controls were 
at a control deflection of 0°. However, if the pitch controls were 
deflected 5° and thus caused an asymmetry in the model, was 

a 

much greater for the model that rolled than for the model that did not 
roll. For any rocket-propelled model test, no exact evaluation of the 
asymmetries can be made since instrument mountings as well as lifting 
surfaces of the model are subject to construction tolerances. For both 
models of this report preflight measurements of wing misalinements showed 
them to be very small and about the same in the X,Y plane as in 
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the X,Z plane . This result is also indicated from the flight measure- 
ments of figure 6. Further investigation will he necessary if the causes 
of this large apparent asymmetry are to he completely understood. 


Damping 

The damping constant could not he reduced directly from the time 
history of a or % since the rather large roll rate caused the 
response to not he confined to the X,Z plane. Two of the methods 
used for reducing data from rolling symmetrical models axe shown in 
references 5 and 6. Both methods involve the making of plots of C]sj 
against Cy or a against p and determining a resultant force coef- 
ficient. In this report the method of reference 5 was used primarily 
hut the method of reference 6 was also used to check agreement. Sample 
plots of the variation of Cjj with Cy are shown in figure 12 for 

model 1 and figure 13 for model 2. These shapes are typical for a small 
value of the space-motion factor K which is described in reference 5 • 

The value of £ for model 1 varied from -0.37 M = 2.^5 to -0.2 
at M = 1.77* The value of £ for model 2 varied from -0.3^ 
at M = 2.25 to -0. 28 at M = I.69. The plots of Cjy against Cy 

from all the pulses were very similar to figures 12 and 13 except for 
the first pulse of model 2. Inspection of the roll-rate trace in 
figure 7 shows that the second pulse rocket fired less than 0.1 second 
after the first pulse rocket finished firing. This model motion was 
somewhat similar to the motion resulting from one long burning pulse 
rocket and results in a very large value of the space-motion factor K. 

In fact, the model had rolled through more than 90° between the start 
of the first pulse rocket to burnout of the second pulse rocket. The 
plot of Cjj against Cy for these pulses, presented in figure 1^, shows 

that the value of K is extremely large. 

The time history of Cr 2 was obtained from the relationship 
Cr2 = ( Cn - CN-t)^ + ^Cy - Cy t ^. A sample time history of as 

obtained from the plot of against Cy in figure 12 is shown 

in figure 15 . The damping constant b as obtained from these time 
histories of Cp2 is presented in figure 16. The damping constant b 
was also determined by using the method of reference 6, and these points 
are also included in figure l6. The excellent agreement between the 
values of b obtained by using the two methods is not surprising even 
with the large value of K for the first pulse of model 2. Close inspec- 
tion of the two methods shows that, when the proper trim points are 
chosen, reference 6 does almost the same thing graphically as reference 5 
does mathematically. 
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The damping derivative + Cm^ shown in figure 17 was determined 

from the damping constant b by using the relationship: 



57*3 Ci^ I Y \ 2V 
mVc J 5 


The values of used in this equation were the root-mean-square 

values of C]y a and Cy^. 


Also included in figure 17 are calculated values of 



+ 0 




obtained by using the method of reference 7* The fact that better agree- 
ment between theory and experiment is obtained for model 2 must be attrib- 
uted to coincidence. The theoretical C m + C„. for model 2 was, of 

xuq 


course, obtained by adding the effect of the forward surface and its 
downwash to the calculated values for model lj therefore, the theory 
must have overestimated the contribution of the forward surface. 


Pitching Moment 

The static stability derivative was obtained from the following 

relationship: 


^ " 57.3qSc 


3 (a) 2 + b 


C 

2j _ 


L_ ^niq 2V 
mVc 


This is the equation for a nonrolling model. However equation (l4) of 
reference 5 for a rolling model reduces to this equation if the magnus 
term is assumed to be negligible and Iy/Iy is assumed to be zero. 

The quantity u> was determined from the period of the time history 


of C R 2 . 


Since the period of 



is one -half that of C R , 


03 = it /Period. The second term of the equation for C which is the 

% 

contribution of the vertical translational degree of freedom was extremely 
small (less than 1 percent) for these heavy models having small lifting 

surfaces far away from the center of gravity. A plot of CL, against 

a 


Mach number for both models is shown in figure 18. The aerodynamic 
center, as obtained from figure 18 and the root -mean- square values of 


Cu a and -Cyp, is presented in figure 1 9- Also included in figure 19 
are the aerodynamic-center positions as calculated from the theory of 
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reference 1. As in the case of the damping, the better agreement 
between theory and experiment in the case of model 2 must be attributed 
to coincidence. 


Drag 

The drag coefficient as a function of Mach number is presented in 
figure 20. Excellent agreement was obtained between as determined 

from the CW Doppler radar unit and C D as determined from the longitudi- 
nal accelerometer in the model. Since the angle of attack was always 
very small, it is essentially zero-lift drag. The angle-of-attack vane 
may have influenced the drag somewhat but its influence is believed to 
be small. In reference 6 when a slightly longer rod was used on the 
flow-direction indicator the drag obtained agreed well with that obt ain ed 
from wind-tunnel tests. Reference 8 also indicates that, when a long 
rod is used, the reduction in nose drag is much less. 


CONCLUDING REMARKS 


An investigation has been conducted at high Reynolds number to see 
whether potential-flow theory will adequately predict the stability deriv- 
atives of a cruciform-missile configuration having an aspect -ratio-5 
unswept tail and a fineness-ratio-20 blunt-nose body. One model also 
had cruciform, aspect-ratio-5. 4, forward surfaces. The two models were 
given pitch disturbances while rolling at 5 to 10 radians per second.. 

There appeared to be regions of separated flow on the bodies of both 
models flight tested above a Mach number of 1.9 even at angles of attack 
of less than 1°. 

The slope of the variation of normal-force coefficient with angle 
of attack agreed very well with that obtained by theory for both models 
but the slope of the variation of side-force coefficient with an gle of 
sideslip was much greater than that of normal-force coefficient or theory 
for both models. The cause of this additional side force is at present 
unknown. 

The damping derivative of the model without forward surfaces was 
greater than that predicted by theory but good agreement between theory 
and experiment was obtained with the model having forward surfaces. 

This result indicates that in this case the theory overestimated the 
damping contribution of the forward surfaces. 

The aerodynamic center of the model with no forward surfaces was 
much farther forward than would be predicted by theory. The fact that 
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better agreement between theory and experiment was obtained for the 
model with forward surfaces also indicates that in this case the theory 
overestimated the destabilizing effect of the forward surfaces. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 11, 1957* 
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Figure 1.- Sketch of models tested. All dimensions are in inches. 




Figure 2.- Sketch of forward surface , antenna, and model wing for models tested. All dimensions 

are in inches. 



CONFIDENTIAL 


Figure 3*- Photographs of models tested. L-92973 
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Figure 5 .- Comparison of measured total pressure and calculated total pressure for model 1 
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Figure 6.- Concluded. 
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(a) Model 1. 

Figure 10.- Variation of side-force coefficient with angle of sideslip. 
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(b) Model 2. 

Figure 10.- Concluded. 
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Figure 11.- Concluded. 
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Figure 13 .- Typical plot of Cjj against 


Cy for model 2. 
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Figure 20.- Variation of drag coefficient with Mach number for models 1 and 2. 
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